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ABSTRACT 

This  thesis  presents  a  survey  of  techniques  which  can  be 
used  to  suppress  sidelobes  in  antenna  array  systems.   Methods 
applicable  to  fixed  arrays  include  direct  null  positioning, 
element  spacing  adjustment,  perturbation  of  complex  pattern 
function  zeroes,  and  pattern  synthesis.   Methods  applicable 
to  signal  processing  arrays  include  performance  optimization, 
which  can  also  be  extended  to  permit  system  self-adaptation 
to  changing  noise  conditions. 
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I.   INTRODUCTION 

This  thesis  is  intended  to  present  a  survey  of  techniques 
for  the  suppression  of  sidelobes  of  antenna  array  patterns. 
The  use  of  multi-element  sensor  structures  in  seismographic , 
acoustic,  and  electromagnetic  systems  greatly  extends  the 
range  of  control  over  the  resulting  signal  transmission  or 
reception  performance.   A  single  element  restricts  process- 
ing capability  to  a  single  signal;  a  system  of  N  elements 
provides  the  opportunity  to  independently  manipulate  N  sets 
of  operating  parameters  in  such  a  way  as  to  improve  various 
suitable  measures  of  system  performance. 

Initially  we  present  the  basic  characteristics  of  arrays 
and  define  two  fundamental  performance  measures;  signal-to- 
noise  ratio,  SNR ,  and  directive  gain,  G.   Methods  of  sidelobe 
suppression  applicable  to  fixed  structures  are  then  discussed, 
followed  by  methods  which  optimize  performance  criteria  in 
processing  systems. 

It  should  be  noted  that  while  the  methods  described 
encompass  a  large  portion  of  current  theory  they  certainly 
do  not  exhaust  the  field.   The  capabilities  of  large-scale 
computer-controlled  signal  processors,  the  multitude  of  use- 
ful array  geometries,  and  the  variety  of  suitable  performance 
criteria  insure  that  this  field  will  remain  a  productive  area 
of  theoretical  and  practical  research. 


II.   BASIC  ARRAY  PRINCIPLES 

An  "array"  may  be  generally  defined  as  an  arbitrary 
collection  of  sensors  arranged  in  an  arbitrary  spatial  pat- 
tern used  to  simultaneously  detect  the  local  states  of  an 
impinging  field  or  superposition  of  fields.   Historically, 
the  concept  of  an  array  has  been  nowhere  as  general  as  this; 
arrays  have  evolved  from  the  extension  to  additional  numbers 
of  elements  of  fixed  geometrical  patterns,  such  as  the  circu- 
lar arrays  of  underwater  sensors,  or  the  linear  and  planar 
arrays  of  radio  antennas.   Theories  of  array  behavior  have 
shown  a  corresponding  evolution  from  specific,  regular 
algebraic  forms,  such  as  the  complex  polynomials  of  Schel- 
kunoff  [13]  or  the  Tchebyscheff  polynomials  of  Dolph  [7] , 
to  the  stochastically  controlled  self-adaptive  processors 
of  Mermoz  [11] ,  Widrow  [17] ,  and  Griffiths  [10] ,  where  the 
very  indef initeness  of  observed  signals  provides  the  means 
of  prescribing  the  appropriate  processing  technique. 

A.   PATTERN  FUNCTIONS 

The  elementary  N-element  linear  array,  Figure  1,  has  a 
total  far-field  radiation  vector 

N-l  -   *a    •    a 

V(G)    =   V    (G)    y~    a      eJnK:clS:Lny    =   V    (G)    S(G) 
o  ^=n      n  o 

n~0 

where  V   is  the  radiation  vector  of  each  element,  a   is  the 
o  n 

complex  excitation  of  the  n-th  element,  S(G)  is  the  space 
factor  of  the  array,  as  defined  by  Ramo  and  Whinnery  [12] , 
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Figure  1.   N-Element  Linear  Array 
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and  k  is  the  wave  number.   The  radiation  intensity  is 

K  =  K  Isl2. 
o 

Thus,  in  systems  where  the  principle  of  superposition  is 

applicable,  the  field  patterns  of  an  array  can  be  arranged 

as  a  product  of  the  element  pattern  times  a  function  of 

angular  direction.   Isl   is  called  the  array  factor. 

If  a   is  a  constant  in  both  amplitude  and  phase,  the 

pattern  has  a  principal  maximum  in  the  X-Z  plane  in  the 

direction  given  by 

nkdsinG  =  0 
G  =  0, 
and  nulls  in  the  directions  given  by 

NkdsinG  =  2nrf,  n<N 

_   .  -1  f2nTA 

G   sin   UkdJ  ' 

Having  chosen  a  specific  array  geometry  and  excitation,  it 

becomes  apparent  that  the  response  of  the  array  in  directions 

other  than  the  central  direction  is  no  longer  under  arbitrary 

control;  an  interfering  signal  may  be  located  between  the 

nulls  resulting  in  degradation  of  the  performance  of  the 

system. 

The  concept  of  the  pattern  function  can  be  extended  to 

two  or  three  dimensions,  using  a  general  coordinate  system 

as  shown  in  Figure  2.   In  the  case  of  regular  planar  arrays, 

the  space  factor  is 

S(G,$)  =^^a   e  jk(md  sinGcos#  +  nd  sinGsintf)) 
m  n  J 


ELEMENT  n 


7 >    Y 


Figure  2.   General  Array  Geometry 
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where  (9,0)  are  the  direction  coordinates  of  beam  orientation, 

and  d   and  d   are  the  element  spacing  in  the  X  and  Y  dir- 
x      y 

ections,  respectively.   The  general  3-dimensional  array  has 

a  space  factor 

rw^  /f>%  -^T    Jkr  cosct 
S(Q,<P)    ~2iane        n    n 

n 
where  coso;  =  sinGsinG  cos((p-(p   )  +  cosGcosG  ,  (Qn>$>   >rn)    are 

the  position  coordinates  of  the  n-th  element  and  cc     is  the 

r  n 

angle  between  the  "look  direction"  (G,0)  and  the  element 
position  vector. 

B.   SIGNAL- TO- NOISE  RATIO  AND  GAIN 

Signal-to-noise  ratio,  SNR ,  is  a  critical  parameter  in 
many  receiving  systems,  such  as  echo-ranging  systems  and 
most  amplitude-modulation  systems.   Similarly,  gain,  or 
directive  gain,  G,  is  a  critical  signal  transmission  param- 
eter.  In  considering  problems  of  antenna  sidelobes,  it  has 
been  suggested  that  off-axis  signals  may  cause  unwanted 
interference.   Conversely,  adjustment  of  the  antenna  pattern 
to  lower  sidelobes  may  degrade  on- axis  performance.   Techni- 
ques of  optimization  to  be  described  will  maximize  SNR  or  G, 
thus  achieving  a  balanced  compromise  between  sidelobe  level 
and  main-beam  degradation.   In  general,  SNR  is  defined  as 

_  output  signal  power  without  noise 
output  noise  power  without  signal 

=  P  /P  . 
&      o 

This  parameter  is  meaningless  without  appropriate  definition 

of  the  processing  system.   A  general  system  is  shown  in 


Figure  3.   The  output  signal  is 

S   =Ts(t+f  )  *  h  (t) 
o  ^—  m      m 

m 

where  the  incident  signal  is  assumed  to  have  the  same  struct- 
ure as  a  function  of  time  at  all  sensors,  T  is  the  time 

m 

delay  of  the  signal  at  the  m-th  element,  h   is  the  impulse 
response  of  the  m-th  element,  and  *  represents  the  convo- 
lution operation.   The  signal  power  is  then 

Ps  =  25  f  Ws  <M»  Hm  U»  <  (6))  e  J«VCn  -  tn)  da)/2  „ 
m  n  /  „ 

where  the  computation  has  been  expressed  in  the  frequency 

domain  with  signal  power  spectral  density  W  (60)  ,  filter 

transfer  function  H  (60)  ,  and  radian  frequency  variable  CO. 

The  noise  output  is 

n   =  5n(t  +  X   )    *  h  (t), 
o   <<=—        m     m v  '  ' 

m 

so  the  resulting  noise  power  is 

P   =yj  A  (C0)H  (C0)H*('j))ej6Y£  "  t    )^/;5/o-n 
o  ^-^—  I      o    m    n       v  m     n  d#y2JT. 

m  n  ' 

-  CD 

An  ultimate  objective  is  to  manipulate  the  gain  and  phase 

shift  of  the  filters  in  a  frequency-independent  manner  so 

the  filter  response  can  be  represented  as 

H(C0)  =  K  ej^mH(C0) 
m       m 

with  H(60)  the  frequency-independent  factor.   Then, 

m  n  -eo 

and 
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i  e 


o    s. 


N 


Figure  3.   General  N-Channel  Linear  Processor 
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m  n  -oo 

Finally,  SNR  may  be  expressed  in  matrix  form  as 

KfGK 


SNR   ~ 


k'yk 


where 


K      = 


K0e^O 
K.e^l 


M-l 


J^M-1 


K*  -  Adjoint    [K] 


G      ~    [g      ]  /    g 

l^mn        ^mn 

Y     =    [y      ]  ,    y 

umn  mn 


The  resulting  formula,  a  ratio  of  two  quadratic  forms,  pro- 
vides unique  properties  for  optimization. 

Directive  gain,  G,  as  defined  by  Silver  [15]  is 

power  density  in  direction  (0  ,(p   ) 
G  "" 


average  power  density 
|jj  f     f   P(9,4>)sin9dGd4? 


o   o 


Since,  for  arrays, 


p(9,$)  =  pQ(e,<J»  |s(q,(J))|2 

where  P  (9,^)  is  the  element  power  pattern  and  S(9,^)  is  the 
space  factor  of  the  array,  the  gain  becomes 
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P(e0,<?>0) 


ift  f  f    P0(e,<J>)    | s <e ,cp)j2  sin  eaea<)5 


¥ht  array  factor  is 

2  ="^,^r  ^  ^jk  (r_coso6^-r_cos<xj 
n  m 


„frt    ,+y.i  2  _vv     jk(r  coso6  -r  cosoc  ] 
S(G,q>)|    "Z^nV     n     n   m     m 


where  the  a   now  represent  the  weighting  factors  and  linear 

filter  of  the  receiving  system, 

a   =  K  e^HCO} 
n    n 

Also,  the  transmitted  power  density  in  direction  (G  , $   )  is 

P(G  .0   )  =  P  (G  ,0  )  |s(G  ,  <b   )  |2 
v  o  ro     o   o  Yo    '    o  *o 


so  the  gain  becomes  (Cheng  [4] ) 
P  (G  ,(h  )  |S(G  ,(b   )|  2 

_    _        O        O    M  O     '  O      M3    ' 

G 


23T    Jt 


^//po(G,1))     |s(G,<f))|     2sinGdGdCf) 
o      o 

^    /r,      /fx  \^^r  j^(r   coscpC     -r   cosoc 

P    (G    ,u)  )>>a    a   eJ     v    n  no      m  mo 

o      o    md  '■*-  *—  n  m 


n  m 


2ir  n 


1       l      f   *    ie\    aaV  VI  jk(r    costX -r    costf.  )     .     rt,rt,A 

4Jt7     /     Po^    '^^^?name  n     m  m' sinGdGd? 


o      o  n  m 

=  K^AK 

t       " 

K   BK 

This  is,  again,  a  ratio  of  quadratic  forms,  with  K  the  vector 
of  element  weights,  and 

2JT  Jt 

i       i_  1       r   f  n    /n    ^\    jk(r   costft  -r   cos#  ) 

B         [bmJ,    bmri     "  ^y   J    PQ(G,(|))eJ         n  n     m  m    . 


mnJ         mn 

o   o 


sinGdGd<£> 


13 


A   =    [a      ]    -   DD 
1   mnJ 


D    = 


jkr   coso6 
eJ      o  00 


jkr..    coso^ 
s        1 
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jkrx„   ^cosM 
iJ      M-l  (M-1)0 


C.   SIDELOBE  SUPPRESSION  PRINCIPLES 

The  most  elementary  technique  of  interference  elimination 
is  to  position  a  null  in  the  interfering  direction.   This 
method  is  restricted  to  situations  in  which  all  other 
noises  are  insignificant  and  where  appropriate  care  can  be 
taken  in  system  alignment  to  allow  for  the  sharpness  of 
nulls.   Positioning  a  null  means  adjusting  the  relative 
element  phases  so  that  'complete  cancelation  occurs  in  the 
null  direction.   Since  maximum  signal  reception  or  radiation 
requires  a  different  phasing  arrangement,  the  technique  leads 
to  sacrifice  in  performance.   The  suppression  of  wide  sectors 
of  noise  generally  requires  more  than  null  positioning.   Here 
it  is  desirable  to  specify  some  acceptable  level  of  gain 
suppression  in  the  sector.   This  result  can  be  achieved  by 
adjustment  of  both  element  gain  and  phase.   An  artificial 
method  of  gain  adjustment  involves  modifying  the  positions 
of  the  elements  in  a  regular  array  structure.   The  complex 
polynomial  technique  of  Schelkunoff  permits  both  direct 
positioning  of  zeroes  corresponding  to  desired  nulls,  and 
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clustering  of  zeroes  to  reduce  sidelobe  levels  over  an 
angular  sector. 

Optimization  processes  permit  use  of  the  total  array 
environment  to  adjust  the  element  filters  to  the  optimal 
values.   Thus,  background  noise,  sectoral  noise,  and  narrow- 
beam  interference  may  all  be  considered.   In  addition,  the 
observation  of  the  noise  environment  at  the  sensors  can 
provide  information  which  can  be  used  in  a  real-time  sense 
to  readjust  the  filters,  thus  permitting  the  system  to  con- 
tinually adapt  itself  to  the  interference  environment. 
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III.   FIXED  ARRAY  SYSTEMS 

Fixed  array  systems  are  those  in  which  the  parameters 
of  the  array,  such  as  element  position,  gain,  and  phase,  are 
determined  and  the  array  is  constructed  according  to  the 
resulting  specifications.   Methods  of  sidelobe  cancelation 
in  this  situation  include  direct  null  positioning,  spacing 
adjustment,  perturbation  of  complex  pattern  function  zeroes, 
and  pattern  synthesis. 

A.   DIRECT  NULL  POSITIONING 

Schelkunoff  [13]  showed  that  pattern  functions  for  linear 
arrays  can  be  represented  by  complex  polynomials  on  the  unit 
circle 

£(Z)  =J~:  a  Zm,  Z  =  e^sin9  -<j>) 
m=0  m 

where  <p   is  a  constant  phase  shift  per  element  and  a   is  the 
complex  weighting  coefficient  for  the  m-th  element.   Zeroes 
of  f (Z)  lie  on  the  unit  circle  and  define  the  positions  of 
nulls  of  the  pattern.   An  M-element  array  will  therefore 
have  M-l  roots  of  f(Z),  determined  by  the  a  ,  and  conversely, 
the  a   can  be  selected  to  independently  position  M-l  nulls. 
The  resulting  design  will  require  individual  element  feeds 
capable  of  both  amplitude  and  phase  control. 

The  concept  of  virtual  couplets  provides  a  null  position- 
ing system  described  by  Davies  [6] .   An  M-element  array  has 
a  pattern  which  can  be  interpreted  by  the  pattern 
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multiplication  principle  as  the  product  of  a  couplet,  or 

two-element  pattern,  and  an  M-l  element  space  factor.   This 

reasoning  can  be  extended  to  its  ultimate  reduction  of  the 
i 

array  to  a  collection  of  couplets.   The  couplets  must  have 
the  same  relative  phasing  at  each  stage.   The  resulting 
structure  for  a  four-element  array  is  shown  in  Figure  4. 
In  this  case,  there  are  three  independently-controlled  phase 
shifts,  which  permit  control  of  three  pattern  nulls. 

B.   SPACING  ADJUSTMENT 

This  method,  due  to  Strait  and  Cheng  [16]  uses  adjust- 
ments in  the  spacing  of  the  linear  array  elements  to  produce 
results  similar  to  amplitude  weighting  of  the  previous 
section.   A  particular  advantage,  however,  is  that  element 
gains  can  now  be  fixed  at  values  satisfying  other  require- 
ments, such  as  fixed-power  matched  drivers.   For  a  linear 
array  of  N  elements  with  N  even  and  cophasal  excitation, 
the  space  factor  is  proportional  to 

N/2 


-    >       I    cos     r—  kdsmQ 

n%      n         LV     2  /  J 


For  a  variation  in  spacing  v, 


N/2 


S      ~     >     I    cos 

v  ^rr       n 

n-l 


fe^Md+Vjsin©] 


2*    Incosr^p)kdsinQ  +    p^zljkvsin^ 

^{(")tCn]f 
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ARRAY 


PHASE  SHIFTERS 


Figure  4.   Null-Positioning  Circuit 
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/2n-3 
cos 


where  y-  kdsinO, 

n    \  2  /  d 

For   c  "^sufficiently    small, 

Since    sin^   ^    Tj> 

It  is  then  apparent  that,  given  a  desired  space  factor  as  a 
cosine  series  with  arbitrary  element  weights,  the  series  S 
can  be  equated  term  by  term  and  the  factors  c   computed  in 

such  a  way  as  to  hold  I   constant.   Such  a  result  would 

J  n 

permit  the  design  of  transmitters  for  optimum,  constant 
operating  power  characteristics,  and  would  also  avoid  the 
complications  of  individual  element  amplification  and  phase 
weighting. 

C.   PERTURBATION  OF  ZEROES 

This  approach  involves  making  small  adjustments  in  the 
locations  of  zeroes  of  the  complex  array  polynomial.   A 
similar  technique  has  been  used  by  Baker  [1] ,  to  modify 
features  of  another  otherwise  suitable  pattern.   Since  the 
space  factor  is 
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N-l 

s(G)  =  J!    (z-zn) 

n=l 


the    array   factor    is 


P(G)    = 


S(G) 


2         N"1 


XT       2e    an  [cosh    aR   -    cos    (Y-bnN 


where   Z  -  e^    ,    ty-  kd(sin9)    -  $ 


and  Z   ~  e   n    J  n  are  the  roots  of  the  polynomial.   Since 

it  is  assumed  that  a   is  not  necessarily  zero,  the  root  at 

n  * 

Z   is  not  necessarily  a  null  of  the  pattern.   The  basic  idea 
is  to  start  with  a  known,  approximately  satisfactory  pattern, 
P  (9) ,  and  to  construct  a  desired  modification  of  the  pat- 
tern, as  shown  in  Figure  5.   Variations  may  now  be  made  in 

a   and  b   to  minimize  the  excess  amplitude  of  the  sidelobes 
n      n  ^ 

in  the  sector  of  interest,  or  to  minimize  the  mean-squared- 
difference  between  the  -desired  pattern  and  the  actual  pat- 
tern.  Normally,  it  would  be  necessary  to  define  the  refer- 
ence pattern  over  the  entire  range  from  9  =  0  to  0  =JI. 
Otherwise,  adjustments  resulting  from  differences  in  a  given 
sector  may  cause  reduction  of  performance  in  the  main  beam. 

D.   PATTERN  SYNTHESIS 

This  procedure  is  also  based  on  the  theory  of  Schelkunoff 
A  desired  space  factor  may  be  defined  by  an  arbitrary 
function,  F(0),  over  the  range  9  =  0  to  9  =3t,  or  its 
equivalent  function  of  y,  f (ty)  ,  where 

^=  kdsin9  -  $ 
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SECTOR  TO  BE  MODIFIED 


Figure  5.   Pattern  Modification 
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If  the  array  has  an  odd  number  of  elements  and  the  coeffici- 
ents equidistant  from  the  ends  of  the  polynomial  are  required 
to  be  conjugates,  the  space  factor  becomes 

N 
S(Y)    =     ^     <A      cosft\Jf+   B      sin   nY) 
n=0 

where  the  unknown  A   and  B   are  related  to  the  polynomial 

n      n  c      x 

coefficients.   Then,  the  desired  space  factor  may  be  expanded 
as  a  Fourier  series  so  that 


f  (Y)  =  Jl   (P   cos  nf     +  <D   sin  nY) 
n=0   n  n 

and  the  coefficients  of  the  two  series  can  be  equated. 
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IV.   LINEAR  PROCESSING  ARRAYS 

Two  types  of  signal  processing  systems  utilize  individual 
filtering,  or  weighting,  of  each  array  element  to  optimize  a 
performance  characteristic  of  the  array.   In  one  method  the 
signal  to  noise  ratio,  SNR,  is  maximized.   Two  general 
approaches  to  SNR  maximization  have  been  considered  by 
Mermoz  [11]  using  N-channel  matched  filters,  and  by  Cheng 
[3,4]  using  the  properties  of  Hermitian  matrices.   In  the 
second  method  the  directive  gain  of  the  array  is  optimized, 
subject  to  constraints  which  may  include  specification  of 
the  desired  location  of  pattern  nulls.   The  objective  of  SNR 
maximization  is  to  obtain  a  set  of  filters  which  result  in 
maximizing  knowledge  of  the  presence  of  a  signal  when  the 
signal  and  noise  are  additive.   The  maximization  techniques 
require  knowledge  of  the  structure  of  noise  in  the  surrounding 
environment;  for  the  matched  filter  approach,  the  noise  must 
be  observable  by  measurement  at  the  array  elements  in  order 
to  determine  its  statistical  properties.   This  has  the  advan- 
tage of  permitting  adjustment  of  the  filters  on  the  basis  of 
the  noise  measurement  if  the  noise  environment  changes, - 
however,  it  also  requires  a  capability  to  process  direct 
noise  information  separately  from  filtered  signal  plus  noise. 
The  matrix  approach  is  equivalent  to  the  matched  filter 
approach  in  receiving  systems; it  permits  constraints  on  null 
locations  to  be  included  in  the  optimization. 
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A.   SNR  OPTIMIZATION 

It  was  shown  that  the  signal-to-noise  ratio  can  be 
represented  by 


SNR  - 


KfGK 
K+YK 


where  K  is  the  vector  of  element  weights  and  G  and  Y  are 
square  matrices  relating  the  relative  element  time  delays 
and  the  filter,  signal,  and  noise  power  spectral  densities. 
Cheng  [3]  has  used  the  resulting  function  to  show  that  if 
G  can  be  represented  by 

t 


G  -  F  F 
o  o 


where 


F   = 
o 


-ikr  cosoi 
e  J      o    oo 

e-jkrlCoso6oi 


e  J   M-l     (K-l) 


then  the  maximum  SNR  is 

SNR     =  F^  Y_1F 
max     o     o 

and  the  maximizing  filter  weights  are 

K*   =  Y-1F  . 
M        o 

The  vector  F   consists  of  the  relative  phase  functions  of 
o  r 

the  signal  at  each  element  of  the  array.   Normalizations 
have  been  assumed  so  that  the  results  are  valid  only  within 
a  constant  factor.   Mermoz  [11]  has  used  the  correlation 
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properties  of  signals  and  noises  at  the  various  elements  to 
arrive  at  the  result  that  the  optimizing  filter,  H  (CO),  is 
given  by 

]>C   (CO)  H  (CO)  =    S  *(0» 
"jj-  mn v    n       m    ' 

where  C    is  the  noise  cross-power  spectral  density  of  the 
mn 

m  and  n-th  elements  and  S   is  the  conjugate  of  the  signal 
spectrum.   It  has  been  assumed  that  the  signal  input  to  the 
filters  has  been  delayed  to  compensate  for  the  relative  delay 
of  the  signals;  this  operation  was  previously  performed  by 
the  vector  F  .   The  result,  then,  for  optimal  performance  is 
that 

a)  The  signal  should  be  delayed  by  a  combination 
of  signal  time  delay  and  noise  correlation; 

b)  The  signal  should  be  attenuated  in  proportion 
to  the  noise  correlation; 

c)  The  result  should  be  linearly  filtered  by  a 
matched  filter,  H  (CO)    -    S*  (CO)  . 

Combining  these  results,  the  complete  filter  system  is  given 

by 

H  (CO)  =  (C-1F  )   S  * 
m  o  m  m 

For  the  M-element  system,  F  provides  the  M  synchronizing 
delays,  C  is  the  MxM  matrix  of  noise  correlations  (in  the 
frequency  domain) ,  and  H  is  the  Mxl  vector  of  channel 
filters.   S   is  lxM.   If  the  noise  environment  is  a  single 
coherent  signal  (jamming)  from  a  certain  direction,  the 
optimization  problem  reduces  to  adding  the  element  outputs 
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in  such  a  way  as  to  maintain  synchronization  of  the  desired 
signal,  while  achieving  maximum  cancelation  of  the  inter- 
fering signal.   A  two-element  monochromatic  situation  as 
shown  in  Figure  6  illustrates  the  application  of  these 
filtering  concepts.   For  the  phase  and  time-delay  reference 
at  element  1,  the  relative  signal  and  noise  delays  at  element 
2  are 

tB  =  i  (d  sin  e  ) 


\   =  §  (d  sin  9) 


If  filters  H,  and  H?  are  assumed  to  consist  of  a  signal- 
synchronizing  delay  followed  by  an  inverse  noise  covariance 


filter,  C   ,  then 


H. 


H. 


e    s 


[■  >] 


-1 
'11 

*-l 

12 


'12 

-1 
'22 


In  this  formulation,  amplitude  normalization  has  been 
assumed.   Since  the  noise  is  coherent  and  monochromatic, 


Cll    C22    1/ 


-1    *-l     iCO  (%     -  %   ) 
C12  "  C21   +  e     N    S 


Accordingly,  with  some  algebra, 


H. 


H, 


1-e  J  v  N     s' 


J^T 


~s  <  l-eJ"  "N 


.m% 


m 
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Figure  6.   2-Element  Filter  System 
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For  signal  and  noise  inputs  A  and  N  at  element  1  and  AeJ   s 
and  NeJ   n  at  element  2,  the  filter  outputs  are 

vx  =  <a+n)  [i  -  e-^<TN  -  V) 

V2  =  (Ae-J"*s  +  Ne-JfflTN)  (e^s)  [l  -  e^N  "  \\  . 

Addition  of  the  filter  outputs  results  in  the  final  output, 

V   =  V,  +  V-  =  2A  -  2A  cosCt>(T  -  T  ) 
o     1     2  N     s 

.=  2A  [1  -  cos C0(Z„  -    Z)]. 

N    s 

Certain  comments  are  appropriate  concerning  the  procedure 
followed  before  evaluating  the  result.   Under  the  assumption 
of  coherent  monochromatic  noise,  the  covariance  matrix,  C, 
is  not  invertible,  since 

det  C  =  0. 
This  implies  that  an  infinite  SNR  is  obtainable,  which  means, 
as  stated  by  Mermoz  [11] ,  that  the  noise  is  completely  can- 
celable.  However,  the  construction  of  phase-shifting 
filters  in  the  desired  structure  can  still  be  carried  out. 
The  procedure  also  assumes  that  relative  attenuation 
between  the  sensors  is  negligible,  which  would  be  approxi- 
mately true  in  practice.   The  resulting  output  signal 
contains  no  noise  component;  therefore  the  process  of 
cancelation  is  valid  regardless  of  the  angles  involved. 
However,  the  signal  amplitude  is  reduced  by  delay  of  the 
noise  relative  to  the  signal.   If  the  noise  and  signal  are 

from  the  same  direction,  LT  =  t  ,  and  the  output  is  zero. 

N     s  v 

In  fact,  due  to  the  broadness  of  the  cosine  function  near 
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t  .  -  L      -   0,  resolution  of  the  signal  and  noise  will  be 
N     s 

difficult  over  a  large  angular  range.   The  resulting  output 
signal,  V  ,  should  be  passed  through  a  matched  filter  in 
the  general  case;  for  monochromatic  signals,  this  would  be 
a  high-Q  resonant  filter. 


B.   GAIN  OPTIMIZATION 

Maximization  of  directive  gain,  G,  is  accomplished 
using  the  same  characteristic  of  quadratic  forms  as  used 
for  SNR.   Thus,  if 

r  =  KfAK 
KTBK 


then 


and 


t  -1 
G     -  D  B   D 

max 


^  =  b"1d' 

where  K   is  the  vector  of  maximizing  filter  weights,  A,  B, 
m 

and  D  are  as  previously  defined.   Since  the  elements  of  B 

jk(r      cos  CL   -r      cosoC     ) 


are  2KS1 


mn 


1       f  f   ^    /r,   ^,\       JMr      cos"    -r      cosoC     )     .    n, ,.,,*, 
■Xjr  J    I    P    (0,9)    eJ     v    n  n      m  m'smGdOd'P 

o   o 

it  is  apparent  that  maximizing  G  is  equivalent  to  maximizing 
SNR  under  the  assumptions  of  monochromatic  operation  with 
isotropic  elements  in  omnidirectional  noise. 

It  may  be  desirable  to  utilize  constraints  in  the  opti- 
mization of  SNR  or  G.  In  the  case  of  transmitting  systems, 
radiation  in  a  certain  direction  may  cause  excessive  inter- 
ference; for  receivers,  a  fixed  noise  source  may  be  more 
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easily  handled  by  simply  requiring  a  null  in  its  i 

rather  than  designing  the  correlation  processes  c 

needed.   Cheng  [3] ,  and  Drane  and  Mcllvenna  [9]  hi 

lined  procedures  to  accomplish  the  constrained  ma: 

based  on  transformation  of  the  A  and  B  matrices  to  an 

abridged  form  which  includes  the  effect  of  the  constraining 

relations. 

C.   ADAPTIVE  SYSTEMS 

The  signal  processing  systems  previously  described  can 
be  designed  for  a  static  mode  of  operation   in  which  the 
filter  weights  are  predetermined.   The  array  may  accordingly 
be  adjusted  for  a  fixed  main  beam  direction  and  fixed  noise 
distribution  or  for  scanning  beams.   The  full  capabilities 
of  the  correlation  principles  used  are  achieved  when  the 
design  is  extended  to  time-varying  situations  in  which  the 
sidelobe  structure  depends  on  the  orientation  of  the  main 
beam  or  the  noise  environment  is  subject  to  random  variation. 
In  such  cases,  the  correlation  measurements  can  be  made  on  a 
real  time  basis,  thereby  keeping  the  system  operation  in 
nearly  optimal  condition.   The  time  varying  characteristics 
are  usually  assumed  to  vary  slowly;  otherwise,  the  statisti- 
cal model  of  the  noise  may  be  over-complicated.   The  main 
problem  in  developing  a  self-adaptive  arrangement  is  con- 
struction of  the  correlation  functions.   One  general  scheme 
for  control  of  the  filters  is  shown  in  Figure  7.   Operation 
of  this  system,  as  described  by  Brennan,  et.  al.  [2] ,  is 
cfefined  by 
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Figure  7.   Adaptive  System  for  SNR  Maximization 
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W  -  G(F  -Z) 
Z  =  Q(V  W  V) 


where  F   is  the  steering  vector  and  Q  is  the  Laplace  trans- 
form of  a  linear  filter,  and  W  is  the  vector  of  weights. 
The  desired  covariance  matrix  is 


m 


.  .    =  E  /v.*v.l. 
ij  \    i      J/' 


N 


But   W  V   -    V 


-     2_     W-V. 
o  ."=r:        i    i 

l    1 


so 


jvVvj  =      E  |v*2wiv.| 


Also, 


MW    ~ 


=      E 


N 


m,  .  w. 
li    l 


N 


A  "^ 


N 


*>        V,      W.  V. 

.%-,      Ill 
1    1 


N        * 

5"     V   ■  w.  v. 
i-fl       N      i    i 


and    since 


m    .    -   E  I  v      v.  I 
m  \    m     i  / 


it  is  clear  that 


32 


E  iv*WTW  =  MW. 

The  weights  are  thus  seen  to  be  controlled  by  the  difference 
between  the  steering  vector  and  a  filtered  function  of  the 
product  of  covariance  and  weights.   In  the  first-order 

case  , 

.  _     1 
Q        T6+   1 

and  in  the  second  order  case 


Q 


WS      +    (tl  +  T2)  S  +  1 
Here  T  represents  the  eigenvalues  of  the  smoothing  process, 
Q,  not  to  be  confused  with  relative  signal  and  noise  delays. 
This  form  of  filter  control  leads  to  the  optimal  weights 
which  maximize  SNR.  '  Widrow  [17]  has  considered  another 
procedure  shown  in  Figure  8  in  which  the  system  adapts  to 
match  a  desired  signal  by  minimizing  the  mean  squared  error. 
The  pilot  signal  can  be  injected  internally,  or  added  to  the 
external  signal  and  noise  field.   Doyle's  [8]  comparison  of 
the  SNR-maximization  and  signal-matching  systems  indicates 
that  the  pilot-signal  technique  requires  high  power  to 
achieve  adequate  "lock-on". 

D.   SPATIAL  SIGNAL  DISTRIBUTIONS 

The  systems  previously  described  have  modeled  signals 
and  noise  as  observable  functions  of  time  for  which  time- 
correlations  can  be  constructed;  the  external,  spatial 
properties  of  the  signals  have  been  assumed  to  be  restricted 
to  time  delays  which  are  explicitly  related  to  spatial 
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Figure  8.   Adaptive  Signal-Matching  System 
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direction.   These  assumptions  are  not  overly  restrictive 
in  normal  applications;  however,  they  may  be  eliminated  by 
application  of  a  more  general  theory  of  space-time  cross- 
correlation  functions.   Such  a  theory  has  been  suggested  by 
Childers  and  Reed  [5] .   Use  of  these  concepts  has  been  con- 
sidered by  Cheng  and  Tseng  [4]  in  the  formulation  of 
signal-to-noise  ratio  for  arrays. 
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